Evidences for bouncing evolution before inflation in cosmological surveys 
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Inflationary cosmology with a preceding nonsingular bounce can lead to changes on the primordial 
density fluctuations. One significant prediction is that the amplitude of the power spectrum may un- 
dergo a jump at a critical scale. In this Letter we propose a phenomenological parametrization of the 
primordial power spectrum in this scenario and confront the jump feature with latest cosmological 
data. Performing a global fitting, we utilize this possibility to derive a novel method for constrain- 
ing bounce parameters via cosmological measurements. Combining the CMB, LSS and SNIa data, 
our result interestingly reveals that a nonsingular bounce, if exists, should be a fast bounce which 
happens at a very high energy scale, as we get an upper limit on the bounce parameters. 



I. INTRODUCTION 

Astronomical observations favor an adiabatic and 
nearly scale-invariant power spectrum of primordial per- 
turbation, which can be realized in an inflation model by 
requiring a cosmological scalar field slowly rolling along 
the plateau of its potential. Based on this, it is usually 
assumed that the primordial power spectrum of cosmo- 
logical perturbation is a fixed power function of the co- 
moving wavenumber k over the range of observable scales 
in the detailed methods of data fitting. However, as advo- 
cated by physics of the very early universe, namely the 
trans-Planckian physics P, bounce cosmologies [H H[ 
and so on, it is possible that the primordial spectrum 
shows local features which cannot be described by the 
usual parameterizationsj^, Q. 

As well-known, an inflation model suffers from the 
problem of initial singularity and thus one cannot use the 
effective field approach to describe the universe at the be- 
ginning of its evolution^. This problem can be circum- 
vented in the framework of bounce cosmology, where the 
big bang singularity is replaced by a nonsingular bounce. 
By virtue of the effective field description, one can real- 
ize a bounce through an effective violation of certain en- 
ergy condition, and thus obtain an inflationary scenario 
with a preceding bounceQ. In this scenario, however, 
the cosmological evolution is usual asymmetric with re- 
spect to the bounce point since of back-reaction of pri- 
mordial perturbation and radiationQ. Therefore, there 
are two plausible scenarios for primordial perturbations 
to exit the Hubble radius, with one being during a phase 
of matter-dominated contraction^, [l(| and the other be- 
ing inflation after the bounce. If the initial perturbation 
originate from a Bunch-Davics vacuum and the contract- 
ing phase connects to an expanding one via a nonsingular 
bounce, the spectra produced in both stages are almost 
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scale-invariant without changes as long as the scales we 
are interested are much larger than the duration of the 
bounce. However, they arc different in their amplitudes, 
and thus their combination yields a nearly scale- invariant 
primordial power spectrum with a jump feature at a crit- 
ical scale. 

In the current Letter, we phenomenologically propose 
a parametrization of primordial power spectrum with 
a jump feature motivated by the bounce inflation sce- 
nario and study the constraints on this feature from cur- 
rent astronomical observations. Compared to the stan- 
dard inflationary paradigm, this model involves two addi- 
tional parameters, ks characterizing the comoving wave- 
length of the universe at the bounce point, and T de- 
scribing the slope of the jump in the power spectrum 
and thus corresponding to the measurement of the dura- 
tion of the bouncing phase. In the literature, the tech- 
nique of Markov Chain Monte Carlo (MCMC) global 
analysis has been widely generalized to constrain non- 
standard inflationary cases, namely non-canonical infla- 
tion models fill. Il2|. trans-Planckian physics flU [Til], and 
so on. It has proven to be very powerful to probe 
the parameter space beyond the standard inflationary 
paradigm. As a consequence, we employ the MCMC 
technique to do a global fitting to constrain parameters of 
bounce inflation and give a comparison to current data. 

The letter is organized as follows. In Section II we 
briefly review the cosmological perturbations generated 
in a bounce model with a matter contraction connecting 
to an inflation, and present a smoothed parametrization 
of the primordial power spectrum of this model with a 
jump feature. In Section III we perform a numerical cal- 
culation and compare the result with the CMB and LSS 
data. We give our numerical results in Section IV before 
concluding with a discussion in Section V. 



II. FORMALISM IN A BOUNCE-INFLATION 
SCENARIO 

We begin with a brief discussion of the cosmological 
evolution of primordial perturbation in the framework of 
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a flat FRW Universe. A standard process of generating 
primordial power spectrum suggests that, cosmological 
fluctuations should initially emerge inside a Hubble ra- 
dius, then leave it in the primordial epoch, and finally 
reenter at late times. In usual, this process can be real- 
ized by stretching the physical wavelength % longer than 
the Hubble radius 1/H as in inflation; while an alter- 
native is to suppress the comoving Hubble radius 1/aH 
shorter than the comoving wavelength 1/k as in matter 
bounce. 

One often uses a gauge-invariant variable the cur- 
vature fluctuation in comoving coordinates, to charac- 
terize the cosmological inhomogeneities. It is associated 
with a canonical variable v — z£, where z = \[2ta with 
e = —H/H 2 . The equation of motion for the Fourier 
mode Vk(v) in the context of standard Einstein gravity 
is given bv[27| 



(k 2 - —)v k = 



(1) 



where the prime denotes the derivative with respect to 
the comoving time rj = J dt/a. 

To perform a specific analysis, one may take the scale 
factor as a(t) = asij^) 1 ^, where the subscript "_B" 
denotes any reference time which will be referred as 
the bouncing point later. For a constant background 
equation-of-state (EoS) w, one obtains 



z 
z 



with v = ±- 
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We assume the cosmological perturbations originate 
from vacuum fluctuations, which suggests 
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when fluctuations are born with \krj\ 1. This is con- 
sistent with the asymptotic solution to Eq.([T]) when the 
last term — is negligible. Therefore, the mechanism of 
generating primordial perturbations requires that the ab- 
solute value of the comoving time is enough large which 
can only be achieved in a contracting or an inflationary 
setup. Another asymptotic solution to Eq.([T]) can be de- 
rived by virtue of the mathematic property of the Bessel 
function, 



c(k)r,-\"\ 



(4) 



in the super- Hubble regime, which implies \krj\ -C 1. 

Now we match the two asymptotic solutions ([3]) and 
(13} at the moment of Hubble crossing \h]\ ~ 1, and thus 
determine the form of Vk on super-Hubble scale, 
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3/2. As a consequence, the primordial fluctuations are 
nearly scale-invariant in both the matter contraction and 
inflationary scenarios. However, e takes the value 3/2 
in the matter contraction but becomes very small dur- 
ing inflation. The amplitude of the primordial spectrum 
would undergo a jump around the scale comparable to 
the bounce scale. A detailed calculation gives the ex- 
pression of the primordial power spectrum for a model of 
bounce inflation as follows 



H 
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k < ks 
k > ks , 



(6) 



where the parameter ks characterizes the comoving 
wavelength of the universe at the bouncing point. Note 
that, when the perturbation passes through the bouncing 
phase, its positive and negative frequency modes could 
be mixed at the transfer surface and thus bring an oscil- 
lating signal around the bounce scale [gj]. However, this 
signal strongly depends on extra parameters introduced 
in specific models and is not quite sensitive to current ob- 
servations. In this Letter we have smoothed this signal 
and study its average effect directly. 

For phenomenological considerations, we would like to 
parameterize the above spectrum by assuming a form 



Pf = Pr, 



P 



inf 



Prr 



1 + tanh[(fc - k B )T] 



(7) 



where Pi n f = and P m is the spectrum of cur- 

vature perturbation before the bounce. Theoretically, 
Pm = -j^T j which is relevant to primordial tensor fluctu- 
ations, should be less than 1% of Pi n f due to the upper 
bound of tcnsor-to-scalar ratio. The power spectrum af- 
ter the bounce Pinf is parameterized as the power low 
format via Pi n f — Ask™ 3-1 , in which A s and n s are the 
amplitude and the spectral index correspondingly, ks 
is a comoving wavenumber relevant to the bounce scale 
and denotes the occurrence moment of the jump feature 
in the power spectrum (in unit of h Mpc -1 ). The pa- 
rameter T characterizes the slope of the jump feature 
(in unit of Mpc) and thus the duration of the bounc- 
ing phase. In the bounce inflation scenario, we have 
ksT ~ HAts- For this description, the phenomenology 
of the scale-invariant primordial power spectrum will be 
recovered when k — ks > and T — > oo. In the following 
we study observational constraints from current data on 
the bounce parameters introduced in Eq.([7]). 



III. CONSTRAINTS FROM OBSERVATIONAL 
DATA 

A. signal in CMB and LSS 



From the definition of the power spectrum, one learns 
that £ ~ fc 3 / 2 !?^! can only be scale- invariant when \v\ — 



To test our model we first consider the current as- 
tronomical observations from the Wilkinson Microwave 
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Anisotropy Probe 7-ycar data (WMAP7)|16| and the 
Sloan Digital Sky Survey (SDSS)H3|. Similar to the 
trans-Planckian physics, the effect brought by a bounce 
is most sensitive to the modes of primordial perturba- 
tion exiting the Hubble radius at earliest time. There- 
fore, the bounce induced jump feature of the primordial 
power spectrum will imprint its effects on the CMB at 
very large length scales, and correspondingly depress the 
CMB anisotropics at large angular scales. Fig. [T] illus- 
trates this impact on the CMB temperature power spec- 
trum. We compare a ACDM power spectrum (black solid 
line) with io b = 0.023, lo c = 0.11, h = 0.71, r = 0.09, 
A s = 2.2 x 10~ 9 , n s = 0.97 and four bounce models 
with different power spectra, namely P m = 5 x 10 -12 , 
T = 2000, k B = 3 X 10~ 4 (red dashed line), P m = 10" 11 , 
T = 4, k B = 8 x 10~ 6 (green dotted line), P m = 10" 10 , 
T = 2000, k B = 8 x 10~ 6 (blue dash dotted line) and 
P m = ID' 11 , T = 2000, k B = 10~ 3 (cyan short dashed 
line) . 



scale of k B . In numerical computation, we illustrate this 
effect by choosing k B = 3 x 10~ 4 and k B = 10~ 3 re- 
spectively. Additionally, the parameter T influences the 
shape of the temperature power spectrum obviously. If T 
is large enough, the jump feature of the power spectrum 
is shown to be very apparent. However, if we choose a 
small value of T, the whole shape of the power spectrum 
becomes smooth relatively and a quite sizable regime of 
this spectrum could be suppressed manifestly. It can be 
seen by the green dotted line with T — 4, which is too 
small to accommodate the WMAP data, and this feature 
could lead to a constraint on T as will be analyzed later. 

The jump feature in bounce inflation can also leave its 
signature on the matter power spectrum on large scales as 
shown in the lower panel of Fig. [TJ Since the signature of 
the bounce effect occurs only at the edge of the observable 
regime of the current data of LSS and accordingly we can 
only obtain a bound instead of an accurate constraint. 
Our global analysis provides a powerful approach towards 
a future detection of bounce cosmologies making use of 
data from more accurate astronomical experiments. 




k (h Mpc ) 



FIG. 1: Signatures of bounce models in the temperature 
power spectrum (top) and the matter power spectrum (bot- 
tom). We compare a ACDM model (black solid line) with 
four bounce models with different parameters as introduced 
in the main context. The red dots are the observational data 
points from CMB (top) and LSS (bottom). 



At large angular scale, the initial condition is given by 
P m rather than P, n /, so the reduction in the amplitude 
of the power spectrum happens beyond the bounce scale 
where k < k B . This effect can be seen by comparing 
the black solid line with the blue dashed line. Moreover, 
the smaller value of P m we choose, the effect of this re- 
duction is more obvious at large scale. The parameter 
k B is determined by the energy scale of the bounce, and 
thus the minimal size of the universe in this scenario. 
From Fig. [TJ we can find that the amplitude of primor- 
dial power spectrum would obtain an increase around the 



B. MCMC Likelihood Analysis 

The MCMC technique is widely applied to give multi- 
dimensional parameter constraints from observational 
data. In the current Letter, we employ MCMC to gen- 
erate a random sample from the posterior distribution 
V{6\x) of a set of parameters 9 given an event x (for us, 
it is the total data set that used), and obtain 



V(6\x) 



V(x\9)V(9) 
JV{x\0)V(0)d0 



(8) 



via Bayes' Theorem, where V{x\9) is the likelihood of 
event x given the model parameters 9 and V{0) is the 
prior probability distribution of obtaining a model pa- 
rameter value 9. The simulated random observations 
from the posterior distribution are of the likelihood sur- 
face, and from this sample, we can estimate the posterior 
distribution of the parameter of interest. 

For our implementation, we use a generalized version 
of the CosmoMC package fig, in which the cosmological 
parameters are: 



(w&,w c , 6 s ,T,7i s , A s ,P m ,k B ,T) 



(9) 



where uib = and ui c = £l c h , in which ilj, and il c 

are the baryon and cold dark matter densities relative to 
the critical density, Q s is the ratio (multiplied by 100) 
of the sound horizon to the angular diameter distance at 
decoupling, and r is the optical depth to re-ionization. 
The remaining parameters are related to the primordial 
power spectrum given by Eq.([7]). For simplicity, we as- 
sume a purely adiabatic spectrum of fluctuations and a 
flat universe with a cosmological constant whose EoS is 
w = — 1 at the initial moment. 
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In our analysis, we have included the WMAP7 tem- 
perature and polarization power spectra with the rou- 
tine for computing the likelihood supplied by the WMAP 
team, the matter power spectrum from the "Luminous 
Red Galaxies" sample from the SDSSI17f as well as the 
SNIa "Union" compilation (307 sample) jl9l] . In calculat- 
ing the likelihood from SNIa we have marginalized over 
the "nuisance parameter" [2^ . Furthermore, we make use 
of the Hubble Space Telescope (HST) measurement of 
the Hubble parameter H = lOO/i km s _1 Mpc _1 by ap- 
plying a Gaussian likelihood function centered around 
h = 0.742 with standard deviation a = 0.038 21 , and 



take the total likelihood to be the products of the sep- 
arate likelihoods of CMB, SNIa and LSS. Alternatively 
defining \ 2 = — 21og£, we get 



2 _ 2 i 2 , 2 

Xtotal ~ XcMB + XsNIa + XhSS 



IV. COMBINED CONSTRAINTS 



(10) 



From numerical computations, we find that the present 
data are far from putting explicit constraint on bounce 
parameters on current obscrvationally accessible scales. 
However, the situation will improve greatly with the ap- 
pearance of more and more accurate CMB experiments 
in the near future. With the current data we can obtain 
upper bounds on P m , ks and a constraint on T. From 
the numerical results, we find that P m < 7.03 x 10 -11 , 
k B < 2.44 x 10" 4 , and 2.63 x 10 2 < T < 7.98 x 10 5 
at 2<t C.L.. The result also shows that a bounce model 
gives a slightly better (lower) x 2 f° r the dataset than 
the ACDM model by 1.4. Accordingly one may conclude 
that a bounce inflation model is quite efficient to explain 
cosmological observations at large scale. 




FIG. 2: A contour plot of H-AtB from combined constraints 
on bounce parameters P m , ks and T from CMB + LSS + 
SNIa data in the frame of bounce inflation. 



Finally, we calculated the combined constraints on the 
inflationary Hubble parameter H and the duration of 
bouncing phase Ats, present the contour plot at la and 



2(7 respectively in Fig. [2j The yellow star represents a 
non- vanishing best-fit point. We also studied the com- 
bined constraints on the inflationary spectral index n s 
and its amplitude A s using global fitting. Although the 
derivation is very limit, we find that a nonsingular bounce 
before inflation leads to a larger red tilt for the spectral 
index and to more power for primordial fluctuations com- 
pared with the ACDM model. 



V. SUMMARY 

Since a nonsingular bounce happened at an extremely 
high energy scale, it can hardly be tested directly by 
experiments. To find evidence for it, we need to know 
its observational consequences. This issue has been dis- 
cussed widely in the literature, and one potential clue 
is to study primordial perturbations. For example, in 
the context of the Pre- Big-Bang scenario [22J and in the 
cyclic/Ekpyrotic model [2j], the resulting cosmological 
perturbation was found to strongly depend on the physics 
at the epoch of thermalization, and thus an uncertainty 
for prediction is involved [24l [2oT |. In this Letter, mo- 
tivated by a combined scenario of matter bounce and 
inflation, we introduce a parametrization of the primor- 
dial power spectrum in a class of bounce cosmologies, 
and introduce two parameters to characterize the rele- 
vant physics. By performing a global analysis, we show 
that the present observations are quite consistent with 
a featureless power spectrum, and consequently put a 
strong constraint on bounce cosmologies, which suggests 
a fast bounce model at very high energy scale. Although 
far from being conclusive, our numerical method can be 
extended to other bounce cosmologies, and thus is rel- 
evant in light of forthcoming astronomical observations, 
such as PLANCK. 

Note added: while this Letter was being finalized, 
a paper appeared pointing out that the scenario of 
bounce inflation could be embedded into loop quantum 
cosmology [26j. and the data analysis in that work shows 
an abnormal enhancement at the first peak of CMB 
anisotropy. We find it more generic to consider the 
bouncing physics in a model-independent way as intro- 
duced in the current Letter. 
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